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We study the Magneto-Electric (ME) effect from the viewpoint of the Berry phase connec- 
tion and quantum adiabatic charge transport . The linear response theory for the electronic 
polarization P c \ can be interpreted in terms of the flux or the fictitious magnetic field related 
to the Berry phase in the generalized momentum space including external parameters. An ap- 
plied magnetic field modifies the spin configuration, induces adiabatic deformation of the Bloch 
wavefunction, and results in the electronic polarization. This provides a new mechanism for the 
gigantic ME effect. For a cyclic change of the applied magnetic field, even a quantized charge 
transport is possible. 
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The interplay between spin and charge degrees of free- 
dom is one of the central issues in the physics of strongly 
correlated electronic systems. From this viewpoint, the 
cross correlation between these two degrees of freedom 
is of particular interests. The Magneto-Electric (ME) ef- 
fect is a phenomenon where a finite magnetization M 
appears under an applied electric field E or an electric 
polarization P is induced by an external magnetic field 
H. 

This effect has been studied for a long term since its 
theoretical prediction 1 and experimental observations 2 
in Cr2C>3. A phenomenological description of this effect 
is given by the free energy containing otijEiHj. The se- 
lection rule for the ME coefficient atj is obtained from 
the group theoretical considerations. In the case of anti- 
ferromagnetically ordered Cr2C>3, the time-reversal sym- 
metry (R) accompanied by the spatial translation (T) 
is broken due its peculiar crystal structure, in contrast 
to usual antiferromagnets. Instead, its AF ordered phase 
is invariant under the product of the time-reversal sym- 
metry (R) and the spatial inversion (J), which does not 
prohibit the finite matrix element otij . Furthermore, due 
to the several other symmetries in antiferromagnetically 
ordered Cr203, there remain only the diagonal matrix 
elements a xx = a yy = a±, a zz = a^. 1 

Microscopic theories of the ME effect in Cr2C>3 based 
on these phenomenological observations have been pro- 
posed. 3-5 According to these theories, the applied electric 
field shifts the oxygen ions O 2- and breaks the crystal- 
lographic equivalence between several sublattices com- 
posed of Cr atoms. For example, the g-factors on two 
Cr sublattices differ from each other due to the crystal 
field produced by the shifted oxygen ions. This difference 
induces a finite net magnetization. 3 A temperature de- 
pendence of a\\ seems to be well explained by this mech- 
anism. 3,5 When we trace these theories the other way 
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around, a microscopic origin for the electric polarization 
induced by the applied magnetic field was attributed to 
the displacement of the ligand anion atoms, namely the 
lattice displacements. 

Instead, we are going to study a Bloch electron's contri- 
bution to the macroscopic electronic polarization (MEP) 
induced by the applied magnetic field. The MEP in di- 
electrics has been formulated only recently, since an elec- 
tron's position operator is ill-defined in the Hilbcrt space 
spanned by the Bloch wavefunctions which obey the peri- 
odic Born- von Karman (BvK) boundary condition. 7 ' 8 In 
1990s, it is gradually recognized that the MEP in ferro- 
electrics is the average value of the total charges trans- 
ferred through the unit area in the system associated 
with the atomic displacement. 8 The atomic displacement 
and/or lattice distortions are induced also by the exter- 
nal pressure, leading to the piezo-electric effects. 9, 10 

The MEP is closely related to the quantized adiabatic 
particle transport proposed by Thouless and Niu in the 
1980s. 11 They considered a ID band insulator with the 
periodic boundary condition. When some external pa- 
rameters are deformed adiabatically, and are put back to 
their initial values, the particle number transferred dur- 
ing this cycle is quantized to be non-zero integer. This 
means that the MEP increases by a finite amount, even 
though the external parameters get back to their initial 
values. Recently this kind of "d.c." responses induced 
by "a.c." impulses, which are generally called Quantum 
Pump, are extensively studied in meso (nano)-scale sys- 
tems, where a cyclic variation of a potential shape in- 
duced by the gate voltage pumps up electrons from one 
end to the other end of the system. 12-14 However we 
should make a clear distinction between the Quantum 
Pump discussed in these mesoscale systems and that of 
the ID band insulator. Namely, the former systems is an 
open system attached to the leads with finite S-matrix 
elements between two channels, while, in the latter sys- 
tem, the wavefunction in one end and that of the other 
do not overlap due to the presence of the gapped re- 
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gion between them. Futhermorc, in the latter system, 
the particle number transferred through this gapped re- 
gion during a cyclic process becomes quantized in the 
thermodynamic limit, 11 ' 15 ' 16 which is not the case with 
the mesoscale systems. 

Both MEP and Quantum Pump are closely related to 
the Berry phase. Berry phase is the quantal phase associ- 
ated with an adiabatic time-evolution of a wavefunction. 
Let an external parameter evolve from \ to A along a 
path T. 17 Then the adiabatically evolved wavefunction 
\^(X(t),t)) is given as follows. 

|*(A,t)> = 

e -/r X ^ XW (* i n 8t (X')|V x ,|* iMt (X')>-<iA^. nst( X( t )))(l) 

where l'I'inst(A)) is an instantaneous eigenstate of the 
corresponding Hamiltonian H(X) , and we take an ini- 
tial wavefunction |\I>(A(tj) = Xi,ti)) as |^i n st(^i))- The 
r.h.s. of eq.(l) is invariant under the U(l) gauge trans- 
formation; |* inst (A)> -> e^)|* inst (A)} where 0(A;) = 0. 
Futhermore, when this parameter evolves along a closed 
loop r cyc , a total phase factor this w.f. acquires, which is 
called U(l) phase holonomy, can be given by the surface- 
integral of the flux over an arbitrary surface S whose 
boundary is the loop r cyc ; 

p -fr cyc (*'nst(A)|V x |*inst(X)>-dA _ i f s B(\)-dS 
e e , 

where the flux B(X) is defined as follows; 

B(X) = -zV x x (* inst (A)|V x |* inst (A)>. (2) 

Since this flux is independent of the gauge choice of 
the instantaneous eigenstate, B(X) is also called fictitious 
magnetic field, while *4(A) = (*i ns t(A)|V^|^I'i ns t(A)) cor- 
responds to its associated vector potential. 

This kind of the flux also appears in the context of the 
MEP and/or the Quantum Pump, since the current is 
represented by the derivative of the phase of the wave- 
functions. Specifically the integrated particle current can 
be viewed as the flux penetrating a plaquette spanned by 
the crystal momentum and the external parameter. Al- 
though this definition of the flux is slightly modified from 
that of cq.(2) (see eqs.(ll,12)), both of these are essen- 
tially the same mathematical objects, which are called 
field strengths in the framework of the gauge theory. 
Futhermore, the point where the energy degeneracy oc- 
curs plays a role of the source of the flux, so called, mag- 
netic monopole (fictitious magnetic charge). Generally 
speaking, a magnetic monopole is an origin of nontriv- 
ial structures of the Fiber bundle associated with phase 
factors of wavefunctions. 18 

Recently the spontaneous (anomalous) Hall effect ob- 
served in ferromagnets are discussed in the context of the 
integer quantum Hall current from the viewpoint of the 
Berry phase. 19 For example, the Hall conductivity o~ xy , is 
given by the total flux penetrating a plaquette spanned 
by x-component and y-component of the crystal mo- 
menta (2D magnetic Brillouin zone). 20 Throughout these 
extensive studies, it turned out that non-coplanar spin 
configurations and/or spin-orbit couplings in magnets 



generate magnetic monopoles and thus nontrivial distri- 
butions of the flux in their crystal momentum spaces, 
which result in their anomalous Hall currents. 

In this paper, we study the MEP induced by a defor- 
mation of background spin configurations in the mag- 
netic materials. Since the spin ordering fields (external 
parameters) in the magnetic systems can be controlled 
by the applied magnetic field, this corresponds to the 
electronic contribution to the ME effect. As discussed 
above, the macroscopic electronic polarization (MEP) 
and the transverse conductivity a xy are closely related. 
Namely, these two physical quantities are the different 
components of the same field in the generalized momen- 
tum space. Through the analogy to the integer quantum 
Hall effect, we can naturally expect that even the quan- 
tized charge transport might be possible during the cyclic 
change of spin ordering fields. In fact, both of them are 
related to the first Chern numbers associated with the 
filled bands. 

The main results obtained in this paper are summa- 
rized as follows. From the viewpoint of the perturbation 
theory, the Magneto-Electric (ME) effect was known to 
be enhanced when the energy denominator is small, i.e., 
when the band gap is small. However we propose, in this 
paper, a new mechanism of the ME effect, where not 
only the band gap reduction but also the U(I) phase 
associated with the magnetic Bloch wavefunction play 
very important roles in determining the magnitude of 
the ME effect. From this viewpoint, we can classify the 
magnetic dielectrics into two categories, i.e., with and 
without the nontrivial structure of the U(l) Fiber bundle 
associated with magnetic Bloch wavefunctions. Former 
category is a good candidate to which our mechanism 
can be applied. In these systems, a nontrivial topolog- 
ical structure is originated from the band crossing lo- 
cated in the parameter space. As for a specific example 
of this new category of magnetic dielectrics, we construct 
a model where such a non-trivial topological structure is 
realized in its generalized momentum space and gigan- 
tic electronic polarizations are induced by the applied 
magnetic field. In fact, their magnitudes amount to the 
order of e/a 2 ~ I [C/m 2 ], where e is the electron charge 
and we take the lattice constant "a" to be about 4A. A 
typical ME dielectric Fes04 shows an electric polariza- 
tion of ~ 1.0 x 10~ 5 [C/m 2 ] under an applied magnetic 
field H = 5 [kOe]. 6 Compared with this, the ME effect 
discussed in this paper has 5 orders of magnitude larger 
(Sec. III). Thus these electronic contribution to the ME 
effects can be hardly neglected compared with that of 
the displacement of the ligand anion atoms and so on, 
although the electronic contribution to the ME effects 
has been believed to be much smaller compared with the 
latter. 

This paper is organized as follows. In section 2, we give 
our general idea for the Magneto-Electric effect (ME) 
based on the geometrical viewpoints and explain how the 
flux (fictitious magnetic field) is introduced in the con- 
text of the spin dependent MEP. In section 3, wc propose 
a model which has a nontrivial topological structure in 
its generalized momentum space and exhibits a gigantic 
ME response. Those who want to get the overview of 
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this paper can skip this section. When a system is ter- 
minated without electrodes, our gigantic ME effect, as 
well as conventional ME effects, suffers from depolariza- 
tion fields. In section. 4, we give an appropriate mean- 
field arguments on this problems in the context of spin 
dependent MEP. Section 5 is devoted to conclusions. 

2. Geometrical viewpoints of the spin dependent 
MEP 

2.1 Macroscopic Electronic Polarization 

In this paper, we study the change of the electronic 
polarization induced by a perturbation SH added to an 
original Hamiltonian H: 
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el,fi 



{X e \^)fj). 
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(* \X cl jV m )(* m \6H\ * ) 
(E m - E ) 2 
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(E m - E ) 2 
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Here we introduced an electronic current operator by 
J C L M = — ie[H, X e \ tli \. As shown in Appendix A, this ex- 
pression for the electronic polarization SP e i can be also 
derived as an integrated electron current over a period 
during which the perturbation SH is introduced. In the 
followings, we use eq.(6) instead of eq.(4). This is be- 
cause the current operator J e i iM is well-defined on a peri- 
odic lattice and thus easy to deal with, while the position 
operator X^^ is not. 

Let the Hamiltonian depend on some external param- 
eters ip, i.e. H(tp). Then, by considering SH in the above 
formalism as H(ip + 50) — H(<p), we obtain a derivative 
of the electronic polarization P c \ with respect to (p as 



dP 



el,f 



-Y 



(* |Jel,,|*m)(*m|^|*0) 

(E m - E ) 2 



c.c. 



2.2 spin dependent MEP and the Flux 

Now we study the spin dependent MEP. In this case, 
we regard (p as the spin ordering field </>i.sp as follows. At 
first we decouple the on-site Coulomb repulsion by using 
the Stratonovich-Hubbard field (f> and replace it by its 
saddle point solution 4>i t sp as, 



Z[h] 



J VC^VCexp - J L($ SP ,C\C,h) 



L(^s. P .,C\C, h) = 



(3) 



4 

i 

Hm.f 



H 



M.F. 



Here A e i,^ is the sum of all electron's position opera- 
tors and V is the volume of the system. (X e i^) ^ is an 
expectation value of X e \^ before the perturbation is in- 
troduced. According to the Kubo formula for the linear 
response, <5P e i,M = ($Pei,n) h+sh is given as 



E 



3,a 



U 

~2 



E< 



(8) 



where SH = i[H,SH] and \^ m ) is an eigenstate of the 
Hamiltonian H with its eigenenergy E m . Here we assume 
a finite energy gap between the ground state |\fo) and 
the excited state \^> m ) (m > 1) and focus on the zero- 
temperature case. When we transfer the time derivative 
from SH to X cl ^; 

(*o |*ci,„ | * m )(* m \[H,6H]\* ) 

= -(* \[H,X el J\* m )(* m \6H\* ), (5) 

we obtain the following expression for SP c i, 

*0|4,,.|*rr»><*m|^|*0> 



where [er„] [y = x, y, z) are the 2x2 Pauli matrices. This 
saddle-point approximation is nothing but the mean field 
theory and we identify this saddle point solution 4>i,sp 
(4^ith the spin ordering field. We can control this spin 
ordering field by changing h, Since the saddle-point so- 
lution is determined as a function of an applied magnetic 
field (h) . We study the MEP induced by an adiabatic de- 
formation of this spin ordering field </>i,sp, which we will 
call ip in the followings. 

When we consider a commensurate magnetic order, 
the wavefunction of an electronic state is given by the 
Slater determinant composed by the magnetic Bloch 
wavefunctions; 



(R h a, (j,, a\\$ n %0 )= e ik - Rl {a, fi, a\n), 



(9) 



where (a, fi, a\n) = (a, fi, a\n(k, 0)) denotes a periodic 
part of the magnetic Bloch wavefunction and Ri,a,/j, 
(6)and a represent the magnetic unit cell, the sublattice, 
orbital and spin index, respectively By using this single- 
particle magnetic Bloch wavefunction, the electric linear 
response given in eq.(7) reduces to a more compact form 
in terms of Thouless-Kohmoto-Nightingale-Nijs (TKNN) 
formula, 20 

dP e i,v s = _ . 1 f 



dip 



-ze 



dk 



<„ 


dHM d t^[n)(m 


OHm.f. (k,cp) 
80 






2 



c.c. 



/ dn 
\~d~K 



dn 
dpi 



c.c. 



where H(k, 0) is the Hamiltonian in the momcn- 
l^rijn representation and the periodic part of the 
Bloch w.f. \n) = \n(k, 0)} is its eigenvector with an 
eigenenergy e„ = e n {k,0). The inner product between 

ft) = 4r(W^)>) and denotes the con - 

traction over sublattice, orbital and spin indice, i.e. 

T, a ^Jw:( n ^' ^1°' M ' a ))(^^ a ' & a \ n $> $)))■ Here 
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(C) 

Fig. 1. (a) A shadowed rectangular area is spanned by Sip and the 
line parallel to fc-axis, i.e., 8<px [— — , -^—\. n, Sip and dS form the 
right hand coordinate. The gray arrows represent the flux defined 
in eq. (11). (b) A cyclic deformation of the external parameters 
ip in ipi — if 2 space is defined by r cyc . (c) Shaded regions are un- 
physical parameter regions (metallic regions) where a band gap 
closes, while, the white area denotes a gapped region. A cross (x) 
denotes the parameter point where the H.V.B. and L.C.B. forms 
doubly degenerate point (magnetic monopolc) . These parameter 
points arc included in the metallic regions. The MEP induced by 
the adiabatic change in the gapped region (denoted by broken 
line) is governed by the doubly degenerate points located in the 
metallic region. 

the abbreviation "C.B.", "V.B." and "M.B.Z." repre- 
sent the conduction bands, the valence bands and the 
magnetic Brillouin zone, respectively. 

When we take ^ as a two-component vector, i.e. 
(p = (ipi,ip2), and consider the one-dimensional system, 
we can introduce a 3-dimensional generalized momentum 
space spanned by the crystal momentum k and these ex- 
ternal parameters pi and p2 ■ In the following, the range 
of the magnetic Brillouin zone (M.B.Z.) is denoted as 
[ — a^'aH' where a m is a lattice constant of the mag- 
netic unit cell. In this generalized momentum space, we 
next introduce the following vector fields for each energy 



band n, 

B n {k,ifi,(p2) = V x A n (k,(pi,tp2), (11) 
A n (k,(pi,f 2 ) = -i(n(k,(p 1 ,(p 2 )\V\n(k,(pi,ip2M 2 ) 

where V = {dk,d Vl ,d V2 ). This vector field B n (k, pi, ip 2 ) 
is independent of the following U(l) gauge transforma- 
tion of the magnetic Bloch wavefunction, 

(a,^i,a\n,k,p 1 ,p 2 ) 

- e < €^ 1 ^)(a,M,a|n,fc,^ 1 ,^ 2 ), (13) 

while A n (k,<fi,<fi2) changes by V£(fc, ipi, p>2)- Accord- 
ingly, we call B n (k, p>i,(p2) and A n (k, ipi, p 2 ) as the flux 
and gauge field, respectively. By using these vector fields, 
eq.(10) can be expressed as a surface integral of the flux 
over the rectangular surface spanned by Sp and the k- 
axis, which is drawn as the shadowed region in Fig. 1(a), 

-trr Y] f dS-B n (k,ip 1 ,if 2 ). (1 

n: V ,r>. ^ 1 a m a m J 

The direction of dS in eq.(14) is taken so that n — Sip — 
dS forms the right-hand coordinate where n represents 
a unit vector along the fc-axis, i.e., n — (1,0,0) as in 
Fig. 1(a). 

Our strategy is first (i) to reveal the distribution of the 
flux for every filled energy band and then (ii) to deter- 
mine the optimal direction, Sip op t, which maximizes the 
magnetoelectric response given in eq.(14). For this pur- 
pose, we need to identify the source and/or sink of the 
flux, which we will discuss in the next subsection. 

2.3 Dirac monopole and Quantum Pump 

In the region where the flux is well-defined, i.e. the 
n-th band is isolated from its neighboring bands, there 
are neither source nor sink associated with the flux B n , 
because V • B n = V • (V x A n ) = 0. However in those 
points where the the n-th band forms a degeneracy with 
its neighboring bands, the flux B n becomes ill-defined. 
Therefore, these degeneracy points might become sources 
and/or sinks of the flux. 

As shown in Appendix B, when the n-th band and 
its neighboring band, e.g.(n+ l)-th band, are degenerate 
at (k,ipi,p 2 ) = (^Di ^i.Di </?2,d)j this degenerate point 
becomes in general source and/or sink for the flux B n 
(and also B n+ i) with charge 2-k. We will call these dou- 
bly degenerate points as magnetic (anti-)monopoles or 
Dirac (anti-)monopolcs from now on. Near this doubly 
degenerate point, the effective Hamiltonian for these two 
bands can be in general expanded as follows , 

[H]2 x2 (k,Vi,<P2) 

= (e„(fc D , <Pi,d, <P2,d) + K ■ c)[l] 

+ ^^K,>„], (15) 

li,v 

where K is (k - fc D , Pi - ¥?i,d, P2 - <^2,d), and [<t„] are 
the 2x2 Pauli matrices. The vector c is not important for 
our purpose, while the determinant of the 3x3 matrix 
V plays an important role in classifying this degeneracy 
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point. According to Appendix B, for the flux of the upper 
band ,i.e. B n+ i(k,tp), this degeneracy point acts as an 
anti-monopole (sink) in the case of detV < 0, while as a 
monopole (source) in the case of detV > 0. This relation 
is reversed for the flux of the lower band B n (k, 0), i.e. 

-27r-sign(dctV r )(5 3 (i?), 
27r-sign(dc«0(5 3 (i?). 



V-B„ 
V • B n+1 



(16) 
(17) 



Higher order degenerate points such as fourth-fold de- 
generate points have also a chance to become sources 
and/or sinks for the flux B n . However it often happens 
that these higher order degenerate points in the 3D pa- 
rameter space can be decomposed into several number of 
doubly degenerate points, namely, an effective Hamilto- 
nian reduces to the direct sum of two by two matrices 
given in eq.(15). Therefore we will concentrate on the 
doubly degenerate points in the following. 

Due to the nature of doubly degenerate points given in 
eqs.(16) and (17), when we fill both of these two bands 
by electrons, B n {k,(p) and B n+ i(k, 0) cancel each other 
in the summation of the r.h.s. of eq.(14) and the de- 
generacy point between the n-th and (n + l)-th band 
becomes almost irrelevant for the MEP. On the other 
hand, when only the lower band is filled by electrons, 
B n (k, <p) contributes to the MEP without being cancelled 
by its counter part B n +i(k,0) and results in a gigan- 
tic linear response, where the doubly degenerate point 
between these two bands plays an essential role in the 
direction and magnitude of the MEP. Those magnetic 
dielectrics which belong to the latter case of the elec- 
tron filling are good candidates where this mechanism 
becomes relevant. Based on a specific model, we will dis- 
cuss in section 3 about the physical consequences of these 
two types of electron fillings. 

When we fix the electron number per site, those 
parameter points at which the highest valence band 
(H.V.B.) and the lowest conduction band (L.C.B.) form 
doubly degenerate points are usually unstable points (un- 
physical points). This is because making a direct band 
gap lowers the energy of a ground state wavefunction in 
general. Then the spin ordering field usually takes those 
parameter regions where these degeneracies are lifted. 
We will argue later that this is indeed the case with 
a specific model given in section 3. What is important 
and nontrivial is that the MEP induced by the adiabatic 
change in these gapped phase is still controlled by the dou- 
bly degenerate points hidden in its neighboring unphysical 
parameter points (See Fig. 1(c)). 

When we deform p by a finite amount and put it back 
to the initial value fit, i.e., making a loop r cyc in the 
fi — f2 plane as in Fig. 1(b), the total change of the 
electronic polarization induced by this process is the total 
flux penetrating the cylinder surface spanned by r cyc and 
fc-axis, i.e., r cyc x ~]' 

A cyc -P c l = 



dP cl _ 
-jr-r ■ d(p, 
r u<P 

eve ' 



This is because the surface integrals of the flux B n over 
the k = plane and that over the k = — plane 
cancel each other due to the periodicity. Therefore, if 
this cylinder has some magnetic monopoles and/or an- 
timonopoles of B n inside, we have a chance to obtain a 
non-zero change in the MEP after this cyclic deformation 
along r cyc . This kind of the "d.c." response induced by 
the "a.c." impulse is called Quantum Pump, 11 ' 13 since fi- 
nite amounts of electrons are transported and/or pumped 
up from the one side to the other side of the system after 
one cycle of this adiabatic deformation. 

By using eq.(16) and eq.(17), we can simplifies eq.(18) 
into more intuitive form. Namely, when both the (n + 
l)-th and n-th bands are filled, a pair of magnetic 
monopole(anti-monopole) for the (n + l)-th band and 
anti-monopole (monopole) for the n-th band cancel each 
other in the summation of the r.h.s. of eq.(18). As a re- 
sult of this cancellation, the degeneracy points between 
filled bands do not contribute to the A cyc P e i in eq.(18). 
Then a finite Ar cyc -P e i is originated only from the degen- 
eracy points where the highest valence band (H. V.B.) and 
the lowest conduction band (L.C.B.) touch. Namely, we 
obtain a following simple form for the quantized particle 
transport: 



A cyc -P c l — ^ X 



N 

E 



sign(detVi), 



(19) 



= — y 



dS-B n . (18) 



n:V.B. 1 c y c 



where the index i represents the doubly degenerate point 
between the H.V.B. and the L.C.B.. N is the total num- 
ber of these degenerate points. The 3x3 matrix V given 
in eq.(15) is defined for each doubly degenerate point and 
Vi denotes that of the i-th point. 

If the path r cyc is contractible into one point without 
confronting a metallic region where the H.V.B. and the 
L.C.B. form a band crossing, the left hand side of eq.(19) 
is trivially zero and so is Ar cyc P c i- Therefore eq.(19) in- 
dicates that we must choose the path L cyc to be at least 
noncontractible, in order to get non-zero Ar c C P C \- 

Finally but not the least, we want to mention about the 
stability of the magnetic (anti-)monopole. As is evident 
from eq.(15), we need to tune only 3 real- valued param- 
eters in order to lift the double degeneracy. Therefore, 
when we add another parameter </? pC rturb and constitute 
a four dimensional parameter space, the double degener- 
acy occurs on a line in this 4D space. As a result of this 
line degeneracy appearing in the 4D parameter space, 
these magnetic (anti)monopoles are stable against small 
perturbations. We will also revisit this point in section 
3, where we will generalize the ID studies into the 2D 
case by introducing weak interchain hopping terms and 
demonstrate that the magnetic (anti)monopole observed 
in the ID studies is in fact stable against this interchain 
hoppings. 

3. Gigantic ME Effect — spin-orbit coupling 
model — 

Based on the concept developed in section 2, we will 
construct a model where the global rotation in a spin 
space induces a gigantic electronic polarization. In this 
model, the spin-orbit coupling plays a crucial role, be- 



6 J. Phys. Soc. Jpn. 



Full Paper 



Online-Journal Subcommittee 



Mi 




Maf+ M x 



external magnetic field 

(a) 




/Map L Jm* / B 
/ A f-M„ /A f-M AF /A f-\ 



i I 



(b) 



Fig. 2. (a) The sublatticc magnetization Maf + Mj_ is for the 
A sublattice, while —Maf + M± is f° r the B sublatticc, where 
A and B sublattice arc defined in Fig. (b). The ferromagnetic 
moment M± can be controlled by an applied magnetic field, (b) 
A modulated square lattice with the staggered magnetization on 
A and B sublattice. Grey circles denote the anion atoms while the 
transition metal atoms are located on the lattice. Here we draw 
the ID chains extending along the y-direction as bold lines, while 
the interchain couplings arc introduced as thin lines. Futhermore, 
the anion atoms on the ID chain are shifted along the ^-direction 
uniformly. 



cause the global rotation in the spin space cannot affects 
its electronic state without spin-orbit interactions. 

Let's consider an antifcrromangct on a modulated 
square lattice, where its sublattice magnetizations 
collincar to a particular direction (Maf) cant and ac- 
quire a ferromagnetic moment Mj_ as shown in Fig. 2. 
We can control the direction and magnitude of this fer- 
romagnetic moment M± by using an external magnetic 
field. The hopping integral between the nearest neighbor 
sites depends on the intermediate anion ions, for exam- 
ple, oxygen 2 ~. In Fig. 2(b), we introduce the uniform 
shifts of the anion atoms along the rr-dircction, since our 
spin structure itself does not break the spatial inversion 
symmetry. 

By adding an orbital index on each site in the tight 
binding model, we study the dielectric property of a 
model composed of the c?-orbitals with the spin-orbit 
interaction and the canted antiferromagnetic order. At 
first, we give the result of the 1-dimensional system ex- 
tending along the y-dircction (bold line in the upper 
panel of Fig. 2(b)). Next we give the 2-dimensional gen- 
eralization by adding weak electron transfers along the 



^-direction. This generalization will be justified, since the 
magnetic (anti)monopoles we will observe in the ID stud- 
ies are actually stable against small perturbations such 
as an interchain hopping term. 

The crystal field term Hq.f. and the kinetic energy 
extending along the y-direction i?K,y are given as follows. 
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• (20) 



Here is the energy of the rf-orbital: fj, — 

d x 2_ y 2,d3 Z 2_ r 2,d X y,d yz ,d zx . The transfer integrals are 
composed of the symmetric parts [t v,s ] and antisym- 
metric parts [t v ' a ]. The antisymmetric transfer integrals 
^3)^4) te result from the anion ions' shifts and their signs 
are reversed when the anion ions slide in the opposite di- 
rection (—x direction). Namely, these transfer integrals 
break the spatial inversion symmetry /. Also t\ is in- 
duced by the displacement of the anion ions. An on-site 
spin-orbit coupling is introduced as follows, 
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(21) 



As we explained in section 2, our spin ordering field , i.e. 
the canted AF magnetic moment, enters into our mean- 
field Hamiltonian in the following way: 



H 



MF 



t E a-cj 



fi.a 



(22) 
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Here ipt is taken as (— )Maf + M± when i belongs to 
the A(B)-sublatticc. The two-component ferromagnetic 
moment M± = (M± iX ,M± tV ) can be controlled by the 
applied magnetic field as we mentioned above. 

3. 1 Odd number of electrons filling 

Based on this model Hamiltonian, we consider a situ- 
ation where odd number of electrons are filled per mag- 
netic unit cell. To make the following analysis specific, 
we fill 3 electrons per magnetic unit cell. In order to dis- 
cuss the dielectric property of this electron filling, we 
calculate a following flux B: 



(23) 



where a flux for each energy band, i.e. B n (M±,k y ), is 
defined in cqs.(ll,12), where V in these equations is 
regarded as ( ai ^ ± , d] $ ± , g§-) in the above equation. 
The crystal momentum k y ranges [— ^f-, ^j-]. The dis- 
tribution of this flux in the M± — k y space is given in 
Fig. 3, where we find a source of the flux, i.e., a mag- 
netic monopole at (Mj_, k y ) — (0, 0, 0). This source corre- 
sponds to the doubly degenerate point which the H.V.B. 
(n = 3) and the L.C.B. (n = 4) form at T point in 
the case of M±_ — (see Fig.4(b)). This band touch- 
ing point is nothing but the Kramers doublet, since the 
collinear antifcrromagnet (M± = 0) is invariant under 
the time-reversal operation R combined with the spatial 
translation by a y : {R\a y }. Namely, following two Bloch 
wavefunctions are degenerated at k y = point: 21 

(a, /i, a\n, M± = 0, k y ) = 

b,/3 

where a and b denote the sublattice index. [<x x ] ab in 
eq.(24) exchanges A and B sublattice and thus represents 
the translation by a y . In the case of our electron-filling 
(3 electrons per magnetic unit cell), a finite magneti- 
zation Mj_ always appears and lifts this degeneracy as 
in Fig.4(a,c). This is because making a band gap low- 
ers the total energy of the ground state wavefunction, 
which we can understand by comparing the energy dis- 
persion at Mj_ = and that of finite M±. (compare 
Fig. 4(b) and (a,c)). Thereby the doubly degenerate point 
we mentioned above turns out to locate at an unstable 
parameter point. However, this degeneracy still plays an 
important role for the dielectric property at finite-M±, 
as we explained in section 2. Namely, even when M±_ is 
deformed within the gapped region (finite- M± region), 
the induced electronic polarization P yje i is determined by 
the flux B(M±, k y ) whose nontrivial distribution is orig- 
inated from this doubly degenerate point at M± = 0. 
Futhermore, when M± is deformed along the path r cyc 
which encloses M± = point clockwise (Fig. 3), an in- 
duced electronic polarization A cyc P yje i amounts to +e 
due to this magnetic monopole. These features do not 
alter drastically as far as we fill odd number of electrons 
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Fig. 3. A distribution of the flux B defined in 
eq.(23). t /U = 0.125, ti/U = 0.1, t 2 /U = 0.025, 
t 3 /U = 0.075, U/U,t 6 /U,t 6 /U = 0.05, X/U = 0.025, 
e x 2_ y2 /U = 1.25, e 3z 2„ r 2/C/ = 1.3, e xy /U = 0.05, 
e V z/U = 0.0, e zx /U = 0.05, M A p || [1,1,1] and 
l-^AF = I- We take the M± x -a,xis and Mj_ y -axis as 
(cos 8 cos (f>, cos 8 sin <f>, — sin 8) and (— sin <f>, cos <p, 0) respectively, 
where Maf = (sin0 cos (j>, sin 8 sin</>, cos 6). The vertical axis of 
this figure is the fc^-axis whose range is taken to be [— —^~] 
and we take a y = 10 for visibility. We normalize the vectors to 
be the same in order to show only their directions. 



per magnetic unit cell. Namely, we always find a mag- 
netic monopole or antimonopole at this unstable point; 
M± = 0, which gives rise to nontrivial dielectric proper- 
ties in its neighboring finite- M± regions. 

3. 2 Even number of electrons filling 

On the contrary, when we fill even number of electrons 
per magnetic unit cell, e.g. 4 electrons per magnetic unit 
cell, the flux which describes the dielectric property, i.e., 
B(M±,k y ) = Y^t=i &n(M±, k y ) has no significant distri- 
bution in the M± — k y space. This is because the degen- 
erate point at (M±,k y ) = (0,0,0) also becomes a sink 
of the flux B4. Namely, the direction of B4 becomes op- 
posite to that of B3 and these two flux strongly set off 
each other. As a result of this strong cancellation, the lin- 
ear response dP y ^ e \/dM± itself becomes very small in the 
finite- Mi regions. Futhermore, the quantum pump does 
not occur even when M± is deformed around an M = 
point. This is because, at around M — 0, there always 
remains a finite band gap between the H.V.B. (n = 4) 
and L.C.B. (n = 5). 
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Fig. 4. Energy dispersions of the lower 6 bands, which are mainly 
composed of t2 g orbitals and constitute the lower Hubbard 
bands. Here we take the same parameter values as those in Fig.3. 
The broken lines denote those of the empty bands, while the 
solid lines represent the filled bands (in the case of 3 electrons 
per magnetic unit cell). 



3.3 2- dimensional generalization 

Let us generalize the above ID studies into the 2D case 
with 3 electrons per magnetic unit cell. We introduce the 
interchain hopping integrals, 
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(25) 



Here we ignore the antisymmetric transfer integrals com- 
ing from the anion ions' shifts, because they do not 
change the conclusion essentially as long as they are not 
so strong. A represents the anisotropy in transfer inte- 
grals between x- and y-directions. Then the y-component 
of electronic polarizations dP y _ c \/dM± in cq.(10) requires 
another integral with respect to the crystal momentum 
k x \ 
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217 n=lJtM ± x[-£-,£-] 



In the followings, we will argue that dP y ^ c i/dM± given 
in eq.(26) will exhibit quite similar behaviors as that of 
the previous ID far as the anisotropy parameter 

A is taken to be small. 

When we constitute a four dimensional parameter 
space spanned by k x ,k y and M±, the doubly degener- 
ate point we have observed in the M± — k y space (see 
Fig.3) is expected to form a degenerate line in this 4D 



parameter space. This is because only 3 real-valued pa- 
rameters are sufficient to lift a two fold degeneracy of an 
hcrmitc Hamiltonian. In the context of our model calcu- 
lations, finite k y , M±^ x and M±, y have already lifted the 
double degeneracy and thus a new axis, i.e. fc^-axis, can- 
not lift this degeneracy. In fact, a simple analysis given 
in Appendix C proves that the double degeneracy occurs 
exactly on the line M± = k y — in this 4D parameter 
space. 

As a result of this line degeneracy, the flux for ev- 
ery k x , i.e. ^2 n —i B n (M±, k x , k v ), and its integral over 
the surface 5M± x [-^-,^~], i-c. dP y ^\{k x )/dM±_ ex- 
hibit a similar behaviour as that of J2 n =iB n (M±,k y ) 
and dP e i/dM± in the previous ID studies. This is be- 
cause, for every k x = constant space (a 3D space spanned 
by k y and M±), we always find the doubly degenerate 
point at (Mj_, k y ) = (0, 0, 0), which becomes a source for 
ELiB n (M±,k x ,k y ). 

After taking the integral of dP Vte \(k x )/dM± over k x , 
we calculated the electronic polarization numerically, 
which is shown in Fig. 5. At the center of the flow diagram 
shown in Fig. 5 (a), we find the vortex, which results from 
the doubly degenerate line we mentioned above. In fact, 
the contour integral of this vector field around an arbi- 
trary loop r cyc enclosing this vortex clockwise produces 
+e/& x , where the factor l/a x comes from the integral 
over k x in eq.(26). 

When it comes to the magnitude of dP Vie \/dM±, 
it becomes larger when the direct energy gap be- 
comes smaller. For example, we show the behavior of 
dP yte \(\M±\, 9) /d9 as a function of 9 for several constant 
| Mx | in Fig.5(b), where 9 is defined as (Mj_ >x , M± iV ) = 
|Mj_|(— cos 0, sin 0). In this figure, we can see that 
dPy te \{\M±\,e)/d6 enhances around 9 = 0.5 - 0.6 in 
the case of \M±\ = 0.2. This enhancement is actually 
accompanied by the reduction of the direct band gap. 

When we stack this 2D system in the z-direction with 
weak interlayer couplings, the magnitude of dP y _ c \/d9 
amounts to 1 [C/m 2 ] in the case of a = 4A: 

~ Ji (28) 
89 "a 2 ' 

Meanwhile the electronic energy gain in the pres- 
ence of the finite ferromagnetic moment Mj_ (see 
Fig. 4) is phenomenologically ascribed to the symmet- 
ric part of an anisotropic exchange interactions such as 
Y^i,n= x ,y,z AJ^Si^Si+i^. 22 Thus, as far as the linear re- 
sponse is concerned, we can roughly estimate the mag- 
netic field required to rotate the spin ordering field as 
follows: 

— ~ (29) 

dh~\Aj\- [ yj 

Accordingly, in combination with eq.(28), the magneto- 
electric coefficient in our model turns out to amount to 
10" 3 [C/m 2 • kOe] even in the case of |AJ| = lOmcV. 
This value is in fact much larger that the typical value 
~ 2 x 10~ 6 [C/m 2 • kOe] of the conventional magneto- 
electric materials. On the other hand, when the applied 
magnetic field is sufficiently large compared with this 
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Fig. 5. (a) The flow of dP y (M ± )/dM ± in the 2D case. Here we 
take A = 0.2. The other parameters are taken to be same as 
those of Fig. 3. Here the magnitudes of the vectors are normalized 
to be same, while the informations on their magnitudes arc partly 
given in fig.(b). (b) dP y ([ M± |, 9)/dd as a function of 9, where 8 
is defined in the text. 



anisotropic exchange interactions 23 {gii&h 3> |AJ|), our 
canted component M± follows the direction of the ap- 
plied field, namely M±_ // h. 

Finally, we want to mention about the relation between 
the dielectric property in our models and the conven- 
tional magneto-electric (ME) effect. In the conventional 
ME materials, breaking the IR symmetry by applied 
magnetic fields induce a finite electric polarization. In 
our case, on the other hand, the electric polarization is 
allowed over the entire region in Fig. 5(a). Namely, our 
model with the anion atoms shifted uniformly already 
breaks both the spatial inversion symmetry / and its 
combination with the time-reversal symmetry IR. How- 
ever we might call this phenomenon ME effect in a wider 
sense, since an applied magnetic field produces an in- 
teresting electric response where a cyclic deformation of 
the background spin configuration pumps up electrons 
toward a particular direction. 

4. depolarization fields - mean-field arguments - 

Before rushing to concluding remarks, we want to 
think over the effect of depolarization fields on our 



new mechanism of gigantic Magneto-Electric (ME) ef- 
fect. Namely, throughout this article, we have been as- 
suming implicitly that our system is embedded into a 
closed circuit, where electronic charges transported from 
one end to the other are always short-circuited. How- 
ever, when the system is terminated without electrodes, 
the electrons accumulating at one boundary cause an 
electro-static potential inside the system, which seems to 
push back valence band electrons into the system. Speak- 
ing more generally, when an electric polarization is not 
short-circuited, polarized charges inside the system al- 
ways induce an electric field in the direction reverse to 
this electric polarization. Thereby one might naturally 
expect that this depolarization field strongly reduces the 
Berry phase contributions of magneto-electric responses 
we have discussed so far. 24 In this section, we will argue 
how the magneto-electric response given in cq. (10,14) 
are suppressed in the presence of depolarization field, by 
taking it into account at the mean-field level. 

A depolarization field induced by polarized charges is 
widely known to depend on a shape of the system and 
is sometimes complicate to figure out for an arbitrary 
shape. 25 Thus, in order to understand as simple as pos- 
sible how the depolarization field affects our novel mech- 
anism of ME effect, we will take our system to have the 
most simplest geometry. Namely our system is transla- 
tionally invariant along y and z directions and have only 
two plane boundaries at x = ±-j. Then we can take it for 
granted that both an electric polarization and the depo- 
larization field induced by polarized charges are always 
parallel to the cc-direction. This makes the following argu- 
ments simple. For example, an electric field penetrating 
the plane x' — x, which we will call as Ed{x), consists of 
two parts. One is proportional to an areal charge density 
within the region x < x' < L/2 and the other is to that 
of -L/2 < x' < x: 



E d (x, <p) = - 



Q + (x,0) Q (x,0) 



(30) 



2e 2e 

In the above equation, we have introduced the areal 
charge density Q ± (x, (p) as follows: 

Q+^(x,0) = 



S Ulx<x'<4 



7 <X r <X 



(*o(¥O|JVH*o(0> #A) 



where N? denotes the density operator at f = (x',y', z') 
and S is the total area of our plane boundary. 

By including the electro-static potential due to this de- 
polarization field, we can generalize our previous mean 
field Hamiltonian Hm.f., namely eq.(8), into the follow- 
ing form: 



N 



Hm.f. (<p) = Hm.f. (<p) + e ^2 E ( X J ' ' X J 



(32) 



where Xj denotes the x-component of a position oper- 
ator of the j-th electron. Correspondingly, the previous 
expression for the variation of the electronic polarization 
w.r.t. (p , i.e. eq.(6), is modified into the following equa- 
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The first term gives the Berry phase contribution of 
magneto-electric effect given in eq.(10,14), while the lat- 
ter term brings about a negative feedback effect on this 
topological contributions, which we will argue in the fol- 
lowings. 

From eq.(30), we can easily relate SE d (x) appearing in 
the latter term of eq.(33) with a total electronic current 
which passes through an arbitrary unit area at x' = x 
while the perturbation 80 is adiabatically introduced: 26 



5E d {x) = -—(q + (x,0+S0)-Q+(x,0)) 
eo V / 



J e \(x) 



(34) 



Namely a total electronic charge Q + (x) + Q~(x) is now 
a conserved quantity in absence of electrodes. This cur- 
rent J e \(x) is always periodic w.r.t. a magnetic primitive 
translation vector; J e i(x) — J c \(x + a m ), since we have 
supposed that the perturbation 80 docs not break the 
periodicity of the magnetic unit cell. 27 However its dis- 
tribution within a magnetic unit cell is not necessarily 
uniform in general, which would make following argu- 
ments cumbersome. Thus, we assume, for sake of sim- 
plicity, that its distribution is almost uniform within the 
magnetic unit cell and replace the r.h.s. of eq.(34) by an 
electron current integrated over the x-coordinate; 27 

1 r L/2 

SE d (x) ^SE d = / J cl (a:) dx. (35) 

e Q L J-L/2 

Since this integrated charge current is nothing but the 
change of an electronic polarization <5P c i 8 (see also ap- 
pendix A) , an electric field due to polarized charges now 
turns out to be proportional only to an electronic polar- 
ization; 



SE d 



SPel. 



eo 



(36) 



Admitting this relation between an electronic polar- 
ization and its resulting depolarization field SE d (x), we 
could rewrite eq. (33) into a following self-consistent 
equation for 8P e i; 28 ' 29 



dS ■ B n (k, ipi,ip2) - X SP e i- 



Here we have introduced the electric susceptibility \ as 
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m^tO 



(VPo|Jcl|*m)(*m|*el| fro) 

(E m - E ) 2 
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X c 
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The first term of the r.h.s. in cq.(37) corresponds to 
the Berry phase contributions to the magneto-electric 



response, while the latter term gives a factor to this 
topological contributions. Namely, the above mean- field 
equation for SP e i gives us a following expression for the 
magneto-electric responses in absence of electrodes: 



(33) dtp 
1 



80 



1+ * 277 n^B 



dS-B n (k,<pi,<p®&) 



The electric susceptibility \ in dielectrics ranges from 5 
to 5000. Therefore, the magneto-electric coefficient we es- 
timated in the previous section (see just below eq.(29)) 
decreases by factor 10 _1 to 10~ 3 in presence of dipo- 
larization fields. Futhermore the electric susceptibility 
usually becomes larger when a direct band gap smaller. 
Thus, in the systems terminated without electrodes, our 
Berry phase mechanism of gigantic ME effect seems to 
inevitably meet a severe reduction due to this depolar- 
ization field. Namely, when one look for those mean fields 
near which the H.V.B. and the L.C.B. form magnetic 
(anti)monopoles, not only the Berry phase contributions, 
i.e. X)n V b &n(k, ipi, f2), but the negative feedback fac- 
tor 1 + x also enhances. However, when the system in ter- 
minated with electodes and electronic polarizations are 
measured through (short-circuited) integrated currents, 
the system is in general free from the depolarization field 
mentioned in this section and thereby our Berry phase 
mechanism of gigantic ME effect does not suffer from its 
negative feedback factor. 

5. Conclusion 

In this article, we have studied the new mechanism 
for a gigantic Magneto-Electric (ME) effect. This is due 
to the geometrical Berry phase of the Bloch electrons, 
which is represented by the flux in the generalized mo- 
mentum space including the external parameters such as 
the direction of the sublattice magnetization. We pro- 
posed a specific model, the multi-band models with the 
spin-orbit interaction, where a magnetic monopole and 
an anti-monopole appear in this generalized momentum 
space and determine the distribution of the flux. These 
(anti)monopoles are attributed to the Kramers doublet, 
where the broken time-reversal symmetry in addition to 
the broken spatial inversion symmetry is essential. This 
model calculation indicates that the ME effect is not sim- 
ply determined by the magnitude of the energy gap, and 
the geometrical structure of the U(l) gauge associated 
with the Bloch wavefunctions also plays an important 
role. A comparison with the conventional ME effect 6 sug- 
gests that the magnitude could be <~ 10 3 times larger 
when the monopole's contribution to the electronic po- 
larization survives. Although we do not have an explicit 
candidate for real systems, it would be very interesting to 
look for the material where this mechanism of gigantic 
ME effect works. The detailed band structure calcula- 
tions will be helpful to reveal the distribution of the flux 
in the generalized momentum space, which are left for 
future studies. 
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Appendix A: the integrated electron current 

In this Appendix, we will argue that the change of the 
electronic polarization SP C \.^ given in eq.(6) is identical 
to a total charge transport J c i lA j. 8 This quantity is defined 
as the total current flowing through the system while the 
perturbation SH is adiabatically introduced from t = 
— oo to t = 0: 



el, n 



J, 



<*(i)|Jel,„|*(t)>dt. (Ai) 

) 

Jv 

Here \^f(t)) is an adiabatically evolved ground state 
wavefunction: 

i^|*(t)) = H + 1 (sHe*"-*** + H.c.) |*(t)). (A2) 



Therefore, in the translationally invariant systems, J i, M 
represents the total charges which pass through a unit 
area normal to the fi =constant plane while t £ [— oo, 0]. 
According to the standard time-dependent perturbation 
theory, \$?(t)) can be expanded by the eigenstate |* m ) 
at t = — oo: 

-iE t | \ " |,T/ \„ I „\„—iE m t 



!*(*)) = l*oK 



where its coefficient a m {t) reads: 

J(E m -E Q +uj-i5)s i 
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5H\V )- 



Then the electronic current density is given up to the 
first order in SH as 



l<*(t)|Jel,,x|*(t)> =teA":< 



UA U ) 



*(<"-*)* + c.c, (A-3) 



E m - E +uj - iS 



+ 



(^m|Jcl, M |^0)<^0|^g|^m) 

E m - E - (w - iS) 



), (A-4) 



where we omitted (*o| ^ei^l^o)) because the system is 
an insulator at £ = — oo. Then by substituting eq.(A-3) 
into eq.(A-l), we obtain J e i i(L( in the following way: 

= * J jci, M M - uA w = °)) + c - c -- ( A ' 5 ) 

% [ CO 10 J 

Here we subtracted 2irS{u)j c \^{uj = 0) from J e i,^, be- 
cause we assume that the system remains an insulator 
from t = — oo to t = 0. When we take the static limit 
(oj — > 0) , J e i l([J given in eq.(A-5) is in fact identical to 
<5 Pel given in eq.(6). 



Appendix B: fictitious magnetic charge 

In this Appendix, we will prove that the doubly de- 
generate point defined in eq.(15) is identical to the ficti- 
tious magnetic charge associated with the flux defined in 
eqs.(ll) and (12). Namely, we will show that the surface- 
integral of the flux over an closed surface enclosing this 
doubly degenerate point becomes quantized to be ±2w. 

First we will consider the surface-integral of 
the flux over an infinitcsimally small box v = 
{(K x ,K y ,K z )\(K x ,K v ,K z ) £ [S X ,5 X ] x [S y ,S v ] x 
[— <5 2 ,<y} enclosing this degenerate point, i.e. 



dS-B n {K). 



(B-l) 



Here we take the vector dS to be directed from inside to 
outside of this box. When we define a new coordinate: 







1 




' K x ' 


K> 




sign(dctV r ) 


■[V]- 


Ky 


K _ 




1 




K z 



[T] 



K x 

Ky 

K z 



the 2x2 reduced Hamiltonian given in eq.(15) can be 
rewritten as follows, 

[H] 2x2 = K' x [a x ] + sign(detV) • K' y [a y ] + K' z [* z ], (B-2) 

where we omitted those terms which are proportional to 
a unit matrix. Then, accompanied with this coordinate 
transformation, a following flux can be introduced, 



B' n (K') = -iV gl x {n,K'\Vg,\n,K'), 



(B-3) 



whose relation to the flux in the old coordinate, i.e. 
B n {K), are given as follows: 



BnA K ) 



detr£[T-V*C(*')> 

V 

B'{K'). 



Here a v M denotes the (v,fi)-th cofactor of the 3x3 matrix 
PI- 

We will show in the following that eq.(B-l) can be 
rewritten into the r.h.s. of eq.(B-5) below. The integral 
in cq.(B-l) consists of the surface- integral over the 6 
rectangulars. When we focus on one of them, e.g. the 
rectangular X+ = {{K x ,K y ,K z )\K x = S x ,(K y ,K z ) £ 
[— 5 y ,5 y ] x [— S Z ,S Z ]}, the integral over this rectangular 
can be transformed as follows; 



I Ix 4 



dK y dK z B n 



T(X+) 



dS' ■ B' n (K'). 



(B-4) 



Here the direction of dS' was also taken so that it pen- 
etrates the parallelogram T(X + ) from inside to outside 
of the box T(v). This relation results from the following 
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J-xy 




' T xz 




n,x 


T 


(IKy X 


T 


dK z )- 


n,y 


T 

z v 




T 




B 'n,z 



mathematics; 

dK v dK z B n 

= dK y dK z (a xx B' nx + d yx B' ny + a zx B' nz ), 



= ( 



= dS' ■ B' n , 

where we used the fact that the parallelogram T(X + ) 
in eq.(B-4) is spanned by \T X y,T yy ,T zy f dK y and 
[TxzjTyzjTzz] 1 dK z and that the 3x3 matrix [T] is pos- 
itive definite by definition. The other surface-integrals 
over remaining 5 faces can be also transformed in a simi- 
lar way. As a result, in the if'-coordinate, cq.(B-l) reads: 



cq.(B-l) = 



dT(v) 



dS'-B' n (K'), 



J J^S' ■ B' n (K% (B-5) 



where S is the sphere enclosing the origin, i.e. K' = 0. 

By using the spherical coordinate; (K' x , K' y , K' z ) — 
K'(sm9cos(j), sin#sin</>, cos 9), the eigenvector of 
cq.(B-2) can be given as follows, 



\n,K>) 



COS 2 
Sin i.gi'P signidctV) 
;m fL e -i0-sign(dctV) 

2 a 
— COS t; 



for E n = K' 
for E n = -K' 



Accordingly, the flux defined in eq.(B-3) can be calcu- 
lated easily, 



B'n(K') 



1 -°-°-)%, 



if' 2 sin (T 5(9 
1 



± 



2K' 2 



sign(detV), 



(B-6) 



where the sign + is for the case where the n-th band is 
the upper band and the — sign is for the lower band. 
Substituting eq.(B-6) into the r.h.s. of eq.(B-5), we get 



f f dS-B n (K) 

J Jdv 

= ±2ir ■ sign(detF), 



(B-7) 

where ± correspond to those in cq.(B-6). To summarize, 
in the case of detV < 0, the flux for the upper band sinks 
into K = 0, while that for the lower band streams from 
this origin with 2ir charge. But in the case of ActV < 0, 
the above relations are reversed. 

Appendix C: Maf collinear to an arbitrary di- 
rection 

In this Appendix, we will argue that a double degen- 
eracy always occurs on the line; M± = k y = in the 
k x —k y —M± space. As we argued in eq.(24), the following 
two Bloch wave functions are energetically degenerated: 

(a, fi, a\n(M± — 0, k)) = 

b,0 



since our 2D system with the collinear AF mean field 
is invariant under the time reversal operation combined 
with the the spatial translation, i.e. {i?|a y }. 

In addition to this, the periodic part of the Bloch func- 
tion at (k x ,k y ) and that of (— k x , k y ) are identical to each 
other, which are also energetically degenerated: 

(a,/j,,a\n(M±,k x ,k y )) = (a,/j,,a\n(M±,—k x ,k y )). (C-2) 

We can easily verify eq.(C-2), when we write 
down our Hamiltonian in the momentum represen- 
tation: E[ ff (*.- a 3±)](o,M,a|6, v ,/J)/g i0iMia /£, 6iVi/ ,- Namely, 
[H(k, Ml)] is composed from the following five parts: 
[H] = [H CF ] + [H K , V ] + [H LS ] + [H MF ] + [H K . X ] , 



[H, 



CFJ 



[H KtV } = 2cos(k y & y )[cr x ] ab [t v ' s ] lx „6 a i3 

+ 2 sin(fc y a !y ) [cr y ] ab [t v ' a }^5 al3 , 

[#LS] = 



8 a b[L\\n V [cr\] a p, 

X—x^y.z 



[^mf] ab ( Uuu [O" X\af3 j 

A 

[Hk,x] = 2A cos(fc x a x ) [er x ] ab [t x ' s ]^ v 6 a0 , 

where a, b denote the sublattice indice and /z, v and a, (3 
are those for the orbital and spin respectively. Since we 
didn't take into account the antisymmetric transfer inte- 
grals along the ^-direction, [H(k, M±)] is clearly invari- 
ant under the sign change of k x . Therefore its eigenvector 
at (k x , k y ) is identical to that of (— k x , k y ). 

Eqs.(C-l,C-2) indicate that the following two Bloch 
functions are energetically degenerated along k y = 
when the ferromagnetic moment vanishes (Mj_ = 0); 

(a, fj,, a\n(Mj_ = 0, k x , k y )) = 

E t^U [-ivv] aft (b, M. PHM ± = 0, k x , -kfi}>3) 

b,0 

Appendix D: Maf collinear to the ^-direction 

In this Appendix, we will argue that when Maf is 
collinear to the z-direction, another double degeneracy 
occurs at X point, i.e. (k x ,k y ) — = kx, ir- 

respective of M±. Namely, a doubly degenerate surface 
appears in the 4D parameter space space spanned by 
k x ,k y and Mj_. 

When M± = with Maf collinear to the z-direction, 
the system is invariant under the following symmetry 
operations; 

E, {a y \& y }, IC 2z , {C 2 x|a y }, 

in addition to {i?|a y }. Here I,a y and C 2a; (z) represent the 
spatial inversion, mirror with respect to the y = plane 
and the spatial rotation by n around the a;(z)-axis respec- 
tively. Then the fc-group at X point is composed of a uni- 
tary subgroup Gx = x {E, IC2 Z } and an anti-unitary 
subgroup {i?|a y } x x {E, IC2 Z }- denotes the spa- 
tial translation by Ri = n(& x + & y ) + m(& x — & y ), where 
n and m take an arbitrary integer. Then the irreducible 
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representation of the unitary subgroup x {E.IC^z] 

becomes e x ' Rl Tc lh ({E, IC2 Z }), where Tc lh denotes the 
representation of the monoclinic point group C\h- The 
character tabic for the monoclinic point group is given 
in Table. I, where we only take the double-valued rep- 





E E lC^z IC2z 


r 3 
r 4 


1 -1 -i i 
1 -1 i -i 



Table Dl. Character table for double- valued representation of 
the monoclinic point group Cih- ICiz = EIC'2z 



resentation since we treat the spinful Hamiltonian with 
spin-orbit interations. E in Tabid represents the 2ir ro- 
tation in the spin-i space whose character is always — 1. 
The dimensionality of these two representation are dou- 
bled due to the presence of the anti-unitary subgroup, 
{i?|a y } x Tg x {E, IC2z}- Namely since the Wigner's 
criterion 30 reads 

£ X X (({R\*y}-uf) 

u£Gx,EGx 

= 2^ e 2 ^-^ + ^)(r Clh (i? 2 ) + r Clh (R*cl)), 

Ri 

= 2^2e^< s y+^(T Clh (E)+T Clh (E)) = OPl) 
Ri 

only 2-dimensional representations are allowed at X 
point irrespective of whether Tc lh is chosen to be T3 
or T 4 . Here \ X m eq.(D-l) denotes the character of the 
irreducible representation of the fc-group at X-point. 

These degeneracies still remain even when we violate 
the time-reversal symmetry {i?|a y } by introducing finite 
Mj_. In the case of M± 7^ 0, the system is still invariant 
under the Tg x {E, {i?/C2 Z |a y }}. Then the degeneracy 
at X point becomes doubled due to its anti-unitary ele- 
ment {i?/C 2z |a y }: 

x x (({Ric 2z \z y }-uf) 

u=T Ai x{E,E} 

_ 2 \^ e 2ik x -(Ri+3 y ) _ _27Y j 
Ri 

where 2N is the total number of the elements of the 
unitary k-group x {E, E}. To summarize, we always 

have a double degeneracy at X point regardless of M± . 
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Here the indice i,j are the site indice (not only magnetic unit 
cell indice but also sublatticc indice arc included), fi and v arc 
orbital indice and a, f3 are those for spin. As for the operators on 
the orbital space, when the basis functions for d-orbital indice 
p, and v were taken as \L 2 = 2(2 + 1), L z = 2), |6, 1), |6, 0)(= 
|6,0)*),|6,-1)(= |6,l)*),|6,-2)(= |6,2}*), we replaced [RW] 
in Eq.(D-2) by 

" 1 " 
1 
10 

10 

1 

However, in the convention of this paper, we take the basis func- 
tions for d-orbital as real functions, which means that [flW] 
must be identity. Accordingly, 5x5 matrices [L] in the spatial 
rotational operators should be taken as in eq.(21). As for the 
operators on the spin space, we take S = -|<r. 
C(n, #) -1 (i) in Eq.(D-3) represents the site which transports 
onto the site i when we apply the spatial rotation C(n,9). In 
contrast to the above operations, the spatial inversion / affects 
only on the site index i because the d-orbitals arc always even 
under /, 
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This ground state wave function ^(t)) is now adiabatically 
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N 
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